DESY 11-123 ISSN 0418-9833 

July 2011 



o 

Search for First Generation Leptoquarks in ep Collisions at 

&: HERA 

o 
m 

Oh; HI Collaboration 

(N 
> 

en 

o 



Abstract 



A search for first generation scalar and vector leptoquarks produced in ep collisions is 
performed by the HI experiment at HERA. The full HI data sample is used in the analysis, 
corresponding to an integrated luminosity of 446 pb _1 . No evidence for the production of 
leptoquarks is observed in final states with a large transverse momentum electron or with 
large missing transverse momentum, and constraints on leptoquark models are derived. For 
leptoquark couplings of electromagnetic strength A = 0.3, first generation leptoquarks with 
masses up to 800 GeV are excluded at 95% confidence level. 
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1 Introduction 



The ep collisions at HERA provide a unique opportunity to search for new particles coupling 
directly to a lepton and a quark. An example of such particles are leptoquarks (LQs), colour 
triplet bosons which appear in a variety of beyond the Standard Model (SM) theories [Q]44]|. 
Particle interactions in the SM conserve lepton flavour. If this property is extended to LQ models 
any such particles produced at HERA would decay exclusively into a quark and a first generation 
lepton, namely an electronj or a neutrino. Searches for such signatures have previously been 
performed at HERA [HH71. A dedicated analysis investigating the production of second and 
third generation leptoquarks has also been performed by the HI Collaboration, where the final 
state contains a muon or the decay products of a tau lepton in combination with a hadronic 
system J8l. 

In this paper a search for leptoquarks coupling exclusively to a quark and a first generation 
lepton is performed using the full e ± p collision data set taken in the years 1994-2007 by the 
HI experiment at HERA. The data were recorded with an electron beam of energy 27.6 GeV, 
which was longitudinally polarised up to 38%, and a proton beam of energy up to 920 GeV, 
corresponding to a centre-of-mass energy y/s of up to 319 GeV. The total integrated luminosity 
of the analysed data is 446 pb _1 , which represents a factor of almost four increase with respect to 
the previously published HI results. The presented results supersede those derived in previous 
searches for first generation leptoquarks by the HI Collaboration 0. 



2 Leptoquark Phenomenology and Standard Model Processes 

2.1 Leptoquark production at HERA 

The phenomenology of LQs at HERA is discussed in detail elsewhere 10. The effective La- 
grangian considered conserves lepton and baryon number, obeys the symmetries of the SM 
gauge groups SU(2)l x U(l)y and SU(3)c and includes both scalar and vector LQs. In the 
framework of the phenomenological Buchmiiller-Ruckl-Wyler (BRW) model [91, LQs are clas- 
sified into 14 types ifTOl with respect to the quantum numbers spin J, weak isospin I and chi- 
rality C (left-handed L, right-handed R). Scalar (J = 0) LQs are denoted as Sf and vector 
(J = 1) LQs are denoted Vf in the following. LQs with identical quantum numbers except 
for weak hypercharge are distinguished using a tilde, for example V R and V Q R . Whereas all 14 
LQs couple to electron-quark pairs, four of the left-handed LQs, namely S$, Sf, V Q L and V-y, 
may also decay to a neutrino-quark pair. In particular, for Sq and V L the branching fraction of 
decays into an electron-quark pair is predicted by the model to be (3 e = r eq /(T oq + Y vq ) = 0.5, 
where r eq (T uq ) denotes the partial width for the LQ decay to an electron (neutrino) and a quark 
q. The branching fraction of decays into a neutrino-quark pair is then given by (3 U = 1 — (3 e . 

Leptoquarks carry both lepton (L) and baryon (B) quantum numbers. The fermion number 
F— L + 3 B is conserved. Leptoquark processes in ep collisions proceed directly via s-channel 



'in this paper the term "electron" is used generically to refer to both electrons and positrons, if not otherwise 
stated. 



resonant LQ production or indirectly via w-channel virtual LQ exchange. A dimensionless 
parameter A defines the coupling at the lepton-quark-LQ vertex. For LQ masses M L q below 
y/s, the s-channel production of F = 2 (F = 0) LQs dominates in e~p (e + p) collisions. For 
LQ masses above \/s, both the s and u-channel, as well as the interference with SM processes, 
are important such that both e~p and e + p collisions have similar sensitivity to all LQ types. 



2.2 Standard Model processes 

The search reported here considers final states where the leptoquark decays into an electron and 
a quark or a neutrino and a quark. Such decays lead to topologies similar to those of deep- 
inelastic scattering (DIS) neutral current (NC) and charged current (CC) interactions at high 
negative four-momentum transfer squared Q 2 . The analysis is therefore performed using event 
selections (see section 14.11) similar to those used in inclusive DIS analyses ifTTI and previous 
LQ searches |[5i 

The SM prediction for both NC and CC DIS processes is obtained using the Monte Carlo 
(MC) event generator DJANGOH (HI, which is based on LEPTO [13] for the hard interac- 
tion and HERACLES Ifl4ll for leptonic single photon emission and virtual electroweak correc- 
tions. LEPTO combines 0(a s ) matrix elements with higher order QCD effects using the colour 
dipole model as implemented in ARIADNE |fT5l . The JETSET program |[i~6l is used to simu- 
late the hadronisation process. Additional SM background contributions from photoproduction 
processes are simulated using the PYTHIA IfTTI event generator, with the GRV-G LO [|T8l pa- 
rameterisation of the photon parton density functions (PDFs). All SM expectations are based on 
the CTEQ6m |fT9l proton parton density function parameterisation, which includes only 12% 
of the HI data analysed in this paper, in addition to 30 pb x of ZEUS data. At high Bjorken x, 
the CTEQ6m parameterisations are dominated by data from fixed target experiments due to the 
limited amount of HERA data included. 

Generated events are passed through a GEANT [20] based simulation of the HI appara- 
tus, which takes into account the running conditions of the data taking. Simulated events are 
reconstructed and analysed using the same program chain as is used for the data. 



3 Experimental Method 

3.1 Data sets and lepton polarisation 

The full HI data sample is made up of 164 pb _1 recorded in e~p collisions and 282 pb _1 in 
e + p collisions, of which 35 pb" 1 were recorded at y/s = 301 GeV. Data collected from 2003 
onwards were taken with a longitudinally polarised lepton beam. As leptoquarks are chiral 
particles, these data are analysed in separate polarisation samples, formed by combining all 
data periods with similar lepton beam polarisation P e = (N R — N L )/(N R + N L ), where N R 
(N L ) is the number of right- (left-) handed leptons in the beam. The average polarisation and 
luminosity of all data sets are detailed in tabled] 
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Table 1 : Centre-of-mass energy yfs, average lepton beam polarisation P e and integrated lumi- 
nosity C of the analysed HI data sets. 

3.2 HI detector 



A detailed description of the HI experiment can be found elsewhere [|2T1 . Only the detector 
components relevant to this analysis are briefly described here. A right-handed Cartesian coor- 
dinate system is used with the origin at the nominal primary ep interaction vertex. The proton 
beam direction defines the positive z axis (forward direction). The polar angle 9 and the trans- 
verse momenta Pt of all particles are defined with respect to this axis. The azimuthal angle <\> 
defines the particle direction in the transverse plane. 



The Liquid Argon (LAr) calorimeter |f22j covers the polar angle range 4° < 9 < 154° with 
full azimuthal acceptance. The energies of electromagnetic showers are measured in the LAr 
calorimeter with a precision of a(E)/E ~ 11% /^E/ GeV © 1% and hadronic energy deposits 
with a(E)/E ~ 50%/^E/GeV © 2%, as determined in test beam measurements EUIS). A 
lead- scintillating fibre calorimeter (SpaCal) ll25l covering the backward region 153° < 9 < 
178° completes the measurement of charged and neutral particles. The central (20° < 9 < 160°) 
and forward (7° < 9 < 25°) inner tracking detectors are used to measure charged particle tra- 
jectories and to reconstruct the interaction vertex. The measured trajectories fitted to the inter- 
action vertex are referred to as tracks in the following. The LAr calorimeter and inner tracking 
detectors are enclosed in a superconducting magnetic coil with a field strength of 1.16 T From 
the curvature of charged particle trajectories in the magnetic field, the central tracking system 
provides transverse momentum measurements with a resolution of op T j Pt = 0.002P r /GeV 
© 0.015. The return yoke of the magnetic coil is the outermost part of the detector and is 
equipped with streamer tubes forming the central muon detector (4° < 9 < 171°). In the very 
forward region of the detector (3° < 9 < 17°) a set of drift chambers detects muons and mea- 
sures their momenta using an iron toroidal magnet. The luminosity is determined from the rate 
of the Bethe-Heitler process ep — > ep'j, measured using a photon detector located close to the 
beam pipe at z = —103 m, in the backward direction. 

The LAr calorimeter provides the main trigger in this analysis [11]. In order to remove 
events induced by cosmic rays and other non-ep background, the event vertex is required to 



2 This device was installed in 1995, replacing a lead-scintillator sandwich calorimeter ETI . 



be reconstructed within ±35 cm in z of the average nominal interaction point. In addition, 
topological filters and timing vetoes are applied. 



3.3 Particle identification and event reconstruction 

The scattered electron is identified as a compact and isolated cluster of energy in the electro- 
magnetic part of the LAr calorimeter with an associated track in the inner tracking detectors. 
The hadronic final state is reconstructed using a particle flow algorithm to combine tracks and 
calorimeter deposits not associated to the scattered electron ||26ll27l . The missing transverse 
momentum P™ 1SS , which may indicate the presence of neutrinos in the final state, is derived 
from all reconstructed particles in the event. 

The kinematic quantities in NC events are determined using the electron method Il28l , 
which uses information exclusively from the scattered electron. In CC events, the kinematic 
quantities are determined exclusively from the hadronic final state ||28l . The leptoquark mass 
-^lq = \/Q 2 lv is reconstructed using the measured kinematics of the scattered electron (had- 
ronic final state) in the analysis of NC (CC) topologies, where y is the inelasticity. 



4 Data Analysis 
4.1 DIS event selections 

Neutral current events are selected by requiring a scattered electron with energy E e > 11 GeV 
and Q 2 > 133 GeV 2 . Additionally, a kinematic cut on the inelasticity 0.1 < y < 0.9 
is employed to remove regions of poor reconstruction, poor resolution, large QED radiative 
effects and background from photoproduction processes 0. Background from neutral hadrons 
or photons misidentified as leptons is suppressed by requiring a charged track to be associated 
to the lepton candidate. Energy-momentum conservation requires that T li (E l — Pj) = 2P°, 
where the sum runs over all reconstructed particles, P z is the momentum along the proton beam 
axis and E® is the electron beam energy. Applying the condition T, i (E i — P l z ) > 35 GeV 
further suppresses the contamination from photoproduction background in which the scattered 
lepton is undetected in the backward direction and a hadron is misidentified as an electron. The 
Y>i(E l — P l z ) requirement also further suppresses the influence of radiative corrections arising 
from initial state bremsstrahlung. 

Charged current events are selected by requiring significant missing transverse momentum, 
pmiss > yi GeV, which is due to the undetected neutrino. To ensure a high trigger efficiency and 
good kinematic resolution, the analysis is further restricted to the region 0.1 < y < 0.85. The 
main SM background is due to photoproduction events, in which the scattered electron escapes 
undetected in the backward direction and transverse momentum is missing due to fluctuations 
in the detector response or undetected particles. This background is suppressed by exploiting 
the correlation between P™ 1SS and the ratio V ap /V p JH of transverse energy flow anti-parallel 
and parallel to the hadronic final state transverse momentum vector Ii29ll . 



4.2 Systematic uncertainties 

The experimental systematic uncertainties included in the analysis of the polarised data taken 
in the years 2003-2007 are described in the following. The systematic uncertainties on the 
1994-2000 data are described in the previous HI publication (51. 

In the NC event samples, a systematic scale uncertainty of 1-3% is assigned to the electro- 
magnetic energy measured in the LAr calorimeter, depending on the ^-coordinate of the impact 
position of the scattered electron. A 0.5% component of this uncertainty is considered as corre- 
lated. In addition, an uncorrected uncertainty on the polar angle measurement of the scattered 
lepton of 2 mrad for 9 e > 120° and 3 mrad elsewhere is also included. 

An uncertainty of 2% is assigned to the scale of the measured hadronic energy for events in 
the CC event samples, of which 1% is considered to be a correlated component. In addition, a 
10% correlated uncertainty is assigned to the amount of energy in the LAr calorimeter attributed 
to noise for events in the CC event samples. 

Other experimental systematic uncertainties are found to have a negligible impact on the 
analysis. The effect of the above systematic uncertainties on the SM expectation is determined 
by varying the experimental quantities by ±1 standard deviation in the MC samples and propa- 
gating these variations through the whole analysis. 

The luminosity measurement has an average uncertainty of 3%. The uncertainty on the po- 
larisation measurement is 3.5% and is found to have a negligible effect on the limit calculations 
performed in section [531 

All data sets are compared to a SM prediction based on the CTEQ6m [19] parameterisation 
of the parton densities inside the proton. The uncertainties of this parameterisation are prop- 
agated through the analysis using the full set of eigenvector PDFs, and the effect is added in 
quadrature to the experimental uncertainties listed above. 



5 Results 



5.1 Mass distributions 

Mass spectra of the four HI data sets taken with a longitudinally polarised lepton beam as 
defined in tableQ]are shown in figured! where both the NC and CC event samples are presented. 
The mass spectra of the complete electron and positron HI data sets are presented in figure[2] A 
good description of the HI data by the SM is observed, where the expectation is dominated by 
DIS processes in all event samples, with small additional contributions from photoproduction. 
Since no evidence for LQ production is observed in any of the NC or CC data samples, the data 
are used to set constraints on LQs coupling to first generation fermions. 



5.2 Statistical Method 

For the limit analysis, the data are studied in bins in the Mlq — y plane, where the NC and CC 
data samples with different lepton beam charge and polarisation are kept as distinct data sets. In 
total, iV b i n = 1408 bins are considered, divided equally between the NC and CC event samples. 
For a given bin i, the predicted LQ signal contribution is denoted s t and the predicted number 
of events in the absence of a LQ signal is denoted 6j. The number of events in the presence 
of a LQ signal is thus s» + bi and is obtained as a function of the LQ mass and coupling by a 
reweighting procedure 0. The limits are determined from a statistical analysis which uses the 
method of fractional event counting, optimised for the presence of systematic uncertainties ll30ll . 
For a given leptoquark type, mass and coupling hypothesis, a test statistic X is constructed as a 
fractional event count of all events: 

AW 

X = ^WiUi, (1) 

*=i 

where the sum runs over all bins and n« is the number of events observed in bin i. The weights 
Wi are chosen such that in the presence of a LQ signal the test statistic X is larger than that 
expected from the SM. In particular, bins with a large and positive s, have weights close to one, 
whereas bins with Sj close to zero have weights close to zero. If s, is negative in a given bin 
due to interference effects, the corresponding bin weight is also negative. This has the desired 
effect that an event deficit in such a bin still leads to a X larger than the SM expectation and 
thus is interpreted correctly as a signal contribution. The presence of systematic uncertainties 
may reduce the sensitivity of a given bin. The weight is therefore defined in such a way as to 
ensure that only bins with both a large signal contribution and small systematic uncertainties 
enter with sizeable weights into the test statistic X. This is achieved by defining the weights as 
solutions of the following set of linear equations [|30l : 



Si = fci 



[ Si + bi)Wi + J2 V S Bw i 



+ k 2 



biWi + ^Vgwj 



(2) 



In this analysis, the constants fci and k 2 are set to one, which is the appropriate choice for testing 
signals with a well defined cross section prediction ||30l . The covariance matrices of all bins in 
the presence (absence) of the LQ signal, V$ B (V^), are calculated as: 

V ij = 2^ a ki °kj and yii = 2^ a ki a kj , (3) 

fc=l fc=l 

where afj 3 (<r^) are the one sigma shifts induced from systematic source k to the number of 
events expected in bin i in the presence (absence) of the LQ signal. The sums run over the N sys 
sources of systematic uncertainty. In the case of negligible systematic uncertainties, equation[3] 
is equivalent to the weight definition used in Q. 

Limits are obtained by performing a frequentist analysis of the test statistic obtained from 
the data, X data . For each signal hypothesis, a large number of MC experiments, typically 
0(1O 5 ), are generated by varying the expected number of events s, + bi within the uncertainties. 



Systematic uncertainties are treated as Gaussian distributions and statistical fluctuations are 
simulated using Poisson statistics. For each MC experiment e, a test statistic X e is calculated. 
A probability p data is calculated as the fraction of MC experiments which have X e < X data . The 
LQ hypothesis is excluded at a given confidence level (CL) if p data < 1— CL. In addition to this 
condition, a power constraint is applied OTTl . The power constraint avoids the exclusion of LQ 
signals beyond the sensitivity of the experiment, which may otherwise occur due to statistical 
fluctuations in the data in the opposite direction to that expected from the LQ hypothesis. A 
probability p la is determined as the fraction of MC experiments with X e < X la . Here, X la 
is the value of the test statistic which corresponds to a la downwards fluctuation of the SM. 
It is determined from a second set of MC experiments, where each experiment e is simulated 
in the absence of a LQ signal, that is by simulating systematic and statistical fluctuations of hi. 
The value X la is determined such that the fraction of MC experiments with X e < X la is equal 
to the single-sided la quantile, 15.9%. LQ models are excluded at 95% CL with the power 
constraint applied, if both p data and p la are below 0.05. 



5.3 Limits 

Exclusion limits are first derived within the phenomenological BRW model [9] described in 
section [27T1 Upper limits on the coupling A obtained at 95% CL are shown as a function of 
the leptoquark mass in figure [3l displayed as groups of scalar and vector LQs for both F = 2 
and F — 0. The presented limits extend beyond those from previous leptoquark and contact 
interaction analyses based on smaller HERA data sets by the HI (HEH and ZEUS d[33]l col- 
laborations. For a coupling of electromagnetic strength A = ^47ra cm = 0.3, LQs produced in 
ep collisions decaying to an electron-quark or a neutrino-quark pair are excluded at 95% CL up 
to leptoquark masses between 277 GeV (V^) and 800 GeV (V L ), depending on the leptoquark 
type. 

Within the framework of the BRW model, the S^, 2 LQ decays exclusively to an electron- 
quark pair, resulting in a branching fraction for decays into charged leptons of j3 e = 1.0, whereas 
the Sq LQ also decays to neutrino-quark, resulting in (3 e = 0.5. The HI limits on Sy 2 and Sq 
presented in this paper are compared to those from other experiments in figure 0] Limits from 
the previous HI publication (51 are also shown. Indirect limits from searches for new physics in 
e + e~ collisions at LEP by the OPAL [341 and L3 [1331 experiments are indicated, as well as the 
limits from D0 II361I37H at the Tevatron and from the CMS 1138113911 and ATLAS EOl experiments 
at the LHC. The limits from hadron colliders are based on searches for LQ pair-production and 
are independent of the coupling A, where the strongest current limit for f3 e = 1.0 (/3 e = 0.5) 
scalar LQs is 384 GeV (340 GeV) as reported by the CMS collaboration. For these leptoquark 
masses, this analysis rules out the S^, 2 and Sq LQs for coupling strengths larger than 0.64 and 
0.14 respectively. The HI limits at high leptoquark mass values are also compared with those 
obtained in a contact interaction analysis iflTI . which is based on single differential NC cross 
sections da/dQ 2 measured using the same data. The additional impact of the CC data can be 
seen in the case of the Sq LQ, where a stronger limit is achieved in this analysis, whereas for 
the 5f/ 2 LQ the two analyses result in a similar limit. 

Signatures similar to those expected from LQ decays also appear in supersymmetric models 
with i?-parity violation [|42l . In such models, the production and direct decay of the u 3 L (d^) 
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squark via a A^ (A' llfc ) coupling is equivalent to the interaction of the S^, 2 (Sq) LQ with a 
lepton-quark pair, and as such the results described in the previous paragraph are also valid 
for these squark types, assuming the direct decay dominates. More general limits on squark 
production with i?-parity violating decays are presented in a dedicated HI publication [|43l . 



Beyond the BRW ansatz, (3 e may be considered as a free parameter and the couplings and 
therefore also the branching ratios to electron-quark and neutrino-quark are not necessarily 
equal. By investigating such a model, mass dependent constraints on the LQ branching ratio 
(3 e can be set for a given value of the electron-quark-LQ coupling A e . Excluded regions in the 
/3 e -Mlq plane for three different coupling strengths are shown for a vector LQ with quantum 
numbers identical to V$ in figure [5ta) and for a scalar LQ with quantum numbers identical 
to Sq in figure |5£b). Similar exclusion limits from the Tevatron (D0 Il36ll37l0 and the LHC 
(CMS ll39l and ATLAS RSI ), which do not depend on A e , are also shown in figure [5l For an 
electron-quark-LQ coupling of electromagnetic strength A e = 0.3 the HI limits extend to high 
leptoquark masses beyond the kinematic limit of resonant LQ production, and for most values 
of /3 e extend considerably beyond the region currently excluded by hadron colliders. 



6 Summary 

A search for first generation scalar and vector leptoquarks is performed using the complete HI 
e ± p data taken at a centre-of-mass energy of up to 319 GeV and corresponding to an integrated 
luminosity of 446 pb _1 . The HI data are well described by the SM prediction and no leptoquark 
signal is observed. Limits are derived on 14 leptoquark types and assuming a coupling strength 
of A = 0.3 leptoquarks are ruled out up to masses of 800 GeV, which is beyond the current 
limits from hadron colliders. 
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Figure 1 : The reconstructed leptoquark mass in the search for first generation leptoquarks in the 
2003-2007 HI data, which was taken with a polarised lepton beam. The left-handed electron 
data (a) and left-handed positron data (b) are shown in the top row; the right-handed electron 
data (c) and right-handed positron data (d) are shown in the bottom row. The luminosity C and 
average longitudinal lepton polarisation P e of each data set is indicated. The NC (solid points) 
and CC (open points) data are compared to the SM predictions (histograms), where the shaded 
bands indicate the total SM uncertainties. 
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Figure 2: The reconstructed leptoquark mass in the search for first generation leptoquarks in the 
full HI electron (a) and positron (b) data. The luminosity C of each data set is indicated. The 
NC (solid points) and CC (open points) data are compared to the SM predictions (histograms), 
where the shaded bands indicate the total SM uncertainties. 
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H1 Search for First Generation F = Scalar Leptoquarks 
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Figure 3: Exclusion limits for the 14 leptoquarks (LQs) described by the Buchmuller, Riickl and 
Wyler (BRW) model. The limits are expressed on the coupling A as a function of leptoquark 
mass for the scalar LQs with F = (a) and F = 2 (b) and the vector LQs with F = (c) and 
F = 2 (d). Domains above the curves are excluded at 95% CL. The parentheses after the LQ 
name indicate the fermion pairs coupling to the LQ, where pairs involving anti-quarks are not 
shown. 
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H1 Search for First Generation Scalar Leptoquarks 
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Figure 4: Exclusion limits on the coupling A as a function of the leptoquark mass for the S^, 2 
(a) and Sq (b) leptoquarks in the framework of the BRW model. The parentheses after the 
LQ name indicate the fermion pairs coupling to the LQ, where pairs involving anti-quarks are 
not shown. Domains above the curves are excluded at 95% CL. Limits from the previous HI 
publication (94-00) are also indicated. For comparison, limits from LEP (OPAL and L3), the 
Tevatron (D0) and the LHC (CMS and ATLAS, \fs = 7 TeV data) are shown for comparison, 
as well as constraints on LQs with masses above 350 GeV from the HI contact interaction (CI) 
analysis. 
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Figure 5: Regions of /3 e -M L Q ruled out by the combination of the NC and CC analyses for (a) 
a vector LQ coupling to e + d (with the quantum numbers of the Vq) and (b) for a scalar LQ 
coupling to e~u (with the quantum numbers of the Sq), where only LQ decays into eq and uq 
are considered. Excluded regions at 95% CL are shown as the coloured areas for three different 
values of the electron-quark- LQ coupling A e . Limits from the Tevatron (D0) and the LHC 
(CMS and ATLAS, y/s = 7 TeV data), which do not depend on A e , are also indicated. 
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